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Abstract: HESS J1626— 490, so far only detected with the H.E.S.S. array of imaging atmospheric Cherenkov telescopes, 
HH ■ could not be unambiguously identified with any source seen at lower energies. 

. Therefore, we analyzed data from an archival XMM-Newton observation, pointed towards HESS Jl 626— 490, to classify 

detected point-like and extended X-ray sources according to their spectral properties. None of the detected X-ray point 
sources fulfills the energetic requirements to be considered as the synchrotron radiation (SR) counterpart to the VHE 
source assuming an Inverse Compton (IC) emission scenario. Furthermore, we did not detect any diffuse X-ray excess 
emission originating from the region around HESS Jl 626— 490 above the Galactic Background. The derived upper limit 



for the total X-ray flux disfavors a purely leptonic emission scenario for HESS J1626— 490. 

To characterize the Interstellar Medium surrounding HESS J1626— 490 we analyzed 12 CO(J=1-0) molecular line data 
from the Nanten Galactic plane survey and H I data from the Southern Galactic Plane Survey (SGPS). We found a good 
morphological match between molecular and atomic gas in the —27 km/s to —18 km/s line-of-sight velocity range and 
HESS J1626— 490. The cloud has a mass of 1.8x 4 M and is located at a mean kinematic distance of d = 1.8 kpc. 
Furthermore, we found a density depression in the H I gas at a similar distance which is spatially consistent with the SNR 
G335.2+00.1. 

Therefore, the most likely origin of the VHE 7-ray emission observed with H.E.S.S. is the hadronic interaction of cos- 
mic rays with a moderately dense molecular cloud, which we detected with Nanten. The application of a detailed 
hadronic model for cosmic ray transport and interaction shows that the cosmic rays could originate from the nearby 
SNRG335.2+00.1. 
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■ 1 Introduction ray point sources in the vicinity of HESS J1626-490 and 

searched for possible diffuse excess emissions above the 
HESS J1626-490 is a VHE 7-ray source of unknown ori- expected Galactic background. Furthermore, we present 
gin, which so far could not be identified with any coun- 12 CO(J=1-0) molecular line survey data taken with the 
terpart at lower wavelengths. This object, with an in- Nanten mm/sub-mm observatory to scan for molecular 
trinsic extension of -5 arcmin (Gaussian FWHM), is lo- clouds, and SGPS H I data to characterize the neutral in- 
cated right on the Galactic plane (R.A.: 16 /l 26 m 04 s , Dec: terstellar medium. 
-49°05 / 13 // ) and was detected by H.E.S.S. with apeak sig- 
nificance of 1.5a |1|. These authors measured a power-law - _ _ T „ , A _ 
s P ectrumwitha P hotonindexof2.2±0.1 s , a£ ±0.2 s , s and 2 XMM-Newton X-ray data analysis and re- 
a flux normalization of (4.9±0.9)x 1(T 12 cm- 2 s" 1 TeV -1 sults 
at 1 TeV at energies between 0.5 TeV and 40 TeV 

In this work we analyzed the data of an archival XMM- We anal y zed the EPIC "P n and MOS data of an archival 

Newton observation to search for an X-ray counterpart to XMM-Newton dataset (ID: 0403280201) pointed towards 

HESS 11626-490. Therefore, we classified the detected X- the direction of HESS J 1626-490. After screening for pe- 
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Figure 1 : Adaptively smoothed EPIC-PN image of diffuse 
X-ray flux in the 3-7 keV energy band with a linear color 
scale (arbitrary units). The instrument FoV is shown as 
a dashed circle (green). The large ellipses (black) denote 
the extraction regions for HESS Jl 626—490 and the back- 
ground test, respectively. The large circle (magenta) to the 
east gives the approximate extension of the X-ray scatter- 
ing halo of 4U 1624-490. Excluded regions are shown as 
crossed-out areas (yellow regions crossed with red lines) 



smoothed image of the diffuse X-ray flux in the 3-10 keV 
band. Point-like sources were removed, based on the 99% 
energy containment radius of the PSF. Also, the read-out 
streak caused by the brightest source in the FoV has been 
excluded. The instrumental and particle-induced back- 
ground components were subtracted using a non-X-ray 
background dataset from the calibration database. Apart 
from some contamination towards the east of the FoV due 
to an X-ray scattering halo caused by the bright binary sys- 
tem 4U 1624-490, the level of diffuse X-ray emission ap- 
pears to be homogeneous and flat throughout the FoV. To 
quantify potential diffuse excess X-ray emission from the 
direction of HESS J 1626— 490 and to compare the signal 
to the expected diffuse Galactic background, we extracted 
a spectrum from an elliptical region compatible with the 
intrinsic VHE 7-ray extent, as well as from a test region in 
the southern part of the XMM-Newton FoV (as indicated in 
Fig.Q]). We found no indication for enhanced X-ray emis- 
sion from HESS J1626— 490 compared to the test region, 
and the spectra from both regions are compatible with typ- 
ical diffuse Galactic X-ray emission expected from such a 
region in the Galactic plane (see, e.g., (4)). We derived 
an upper limit for extended X-ray emission connected to 
HESS J1626-490 of F x , e xcess < 4.9 xlO -12 erg cm" 2 
s _1 (1-10 keV), assuming a powerlaw spectral shape with 
index —2. 



riods of enhanced background activity the net exposures for 
pn and MOS are 4.9 ks and 13.2 ks, respectively. 

Using the standard maximum likelihood technique of the 
Science Analysis Software (SAS) we searched for X- 
ray point sources within the 4a VHE 7-ray contours of 
HESS J1626— 490. Twelve such sources were detected in 
this region above a limiting flux of ~2xl0 -14 erg cm -2 
s _1 . We characterized these sources according to their X- 
ray spectral properties as well as their near-infrared coun- 
terpart from the 2MASS all-sky catalog Q. Based on this 
classification scheme (for more details see 0) we calcu- 
lated the unabsorbed fluxes for all detected point sources, 
assuming different spectral shapes and foreground hydro- 
gen column densities, depending on the source class. By 
far the brightest detected X-ray source features a soft ther- 
mal X-ray spectrum with an unabsorbed flux of 7.9 x 10 -12 
erg cm -2 s _1 (0.5 - 10 keV), and can be identified 
with a triplet system of coronally active main sequence 
stars. All other X-ray sources feature unabsorbed fluxes 
<3.5xl0- 13 erg cm" 2 s" 1 (0.5-10keV). 

We also searched for a potential extended and diffuse X- 
ray component related to HESS J1626— 490. O approx- 
imated the morphology of HESS J1626-490 by fitting a 
2-D Gaussian to the VHE 7-ray excess map. The re- 
sulting intrinsic (with the effects of the instrument point- 
spread-function removed) major and minor axes are 0.1° 
and 0.07°, respectively, with a position angle of 3° west 
to north. From the same region we searched for diffuse 
and extended X-ray emission above the astrophysical and 
instrumental background components. Figure Q] shows a 



3 Nanten 12 CO(J=1-0) and SGPS H I data 

To search for molecular clouds, spatially and morpho- 
logically coincident with HESS J1626— 490, we analyzed 
12 CO(J=1-0) molecular line observations performed by the 

4 m mm/sub-mm Nanten telescope, located at Las Cam- 
panas Observatory, Chile 0. For this work the local stan- 
dard of rest velocity (^lsr) range —240 to +100 kms -1 
was searched. 

Figure [2] shows the Nanten 12 CO(J=1-0) image integrated 
over the ^lsr range —31 to —18 kms -1 . In this inter- 
val we found a 12 CO feature partially overlapping with the 
VHE emission. According to the Galactic rotational model 
of (6), this f lsr range corresponds to a kinematic distance 
of 1.5 to 2.2 kpc. Using the relation between the hydro- 
gen column density (Nh) and the 12 CO(J=1-0) intensity 
W( 12 CO), N H = 1.5 x 10 20 [W( 12 CO)/Kkm/s] Q, we es- 
timate the total mass of this cloud at 1.8xlO 4 M for d = 
1.8 kpc within an elliptical region centered at / = 334.78, b 
= 0.00 with dimensions 0.26x0.30 deg. The corresponding 
average density is 2.1 x 10 2 cm -3 . 

In addition to the Nanten 12 CO molecular line data which 
trace the densest regions of the interstellar gas, we also in- 
vestigated the environment of HESS Jl 626—490 using H I 
data from the Southern Galactic Plane Survey (SGPS). Fig- 
ure [3] shows the ^lsr profile for this region integrated over 
the Galactic latitude range —0.11 to 0.24 deg, which is also 
the extent of the VHE 7-ray signal of HESS J1626-490. 
This image shows a region of increased gas density in spa- 
tial coincidence with the CO molecular cloud, most pro- 
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Figure 2: Nanten 12 CO(J=1-0) image of the region around 
HESSJ1626-490 (linear scale in Kkms" 1 ) integrated 
over the ^lsr range —31 to —18 kms -1 . Overlaid are the 
contours of the VHE emission (blue) and of the H I cloud 
discussed in Sect. [3] (green). The dashed circle (light blue) 
denotes the position and extension of SNR G335.2+00.1. 

nounced in the —31 to —23 kms -1 velocity range. The 
contours extracted from an H I image integrated over this 
velocity range are shown in Fig. [2] We estimated the 
mass of the dense HI region coinciding with the VHE 
and 12 CO features using the signal within the same el- 
liptical region as for the CO cloud (Sect. 0. Using the 
relation between HI intensity and column density from 
|8] (X=1.8xl0 18 cm _2 K _1 km/s) we estimated the mass 
of the cloud as 4.9xlO 3 M with an average density of 
60.1 cm -3 . Furthermore, a local HI density depression is 
seen in the center of the image at an angular separation of 
~2V from the CO molecular cloud and HESS J1626-490. 
It is striking that this feature is consistent in position as well 
as in angular extension with the SNR G335.2+00.1 ll9l[T0l. 

4 Discussion 

In this section we discuss potential emission scenar- 
ios giving rise to the VHE 7-ray signal detected with 
H.E.S.S., particularly in the context of the presented multi- 
wavelength data. 

In the case of a leptonic scenario where low-energy photons 
are up-scattered by relativistic electrons via the IC process, 
X-ray emission is expected to accompany the VHE 7-ray 
signal arising from synchrotron cooling of the same pop- 
ulation of high-energy electrons ifTTl . Even from a first 
glance at the flux levels in X-rays (~10 -13 erg cm -2 s _1 



for point sources and ^10 



12 erg cm -2 s _1 for diffuse 
emission) and at VHE 7-rays (~10 -11 erg cm -2 s _1 ), a 
purely leptonic emission scenario seems unlikely. Follow- 
ing lim and assuming a typical Galactic magnetic field 
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Figure 3: H I SGPS Galactic longitude-velocity plot (linear 
scale in K deg) integrated over the Galactic latitude range 
—0.11 to 0.24 deg. The position and extension of the 12 CO 
cloud is denoted by a dashed circle (black). The double 
arrow (green) shows the Galactic latitude and extension of 
SNRG335.2+00.1. 



strength of 3x10 6 G, and a photon index of Y = 2.2 
(as in the VHE domain), we estimate an integrated X-ray 
source flux of-LlxlO -11 erg cm " 2 s" 1 (0.5-10.0 keV). 
This flux is a factor of ~25 more than what we measured 
from all but one point source which makes an identifica- 
tion of HESS J 1626— 490 with any of these sources un- 
likely. As the only sufficiently bright X-ray point source 
in the FoV can be identified with a triplet system of ac- 
tive main sequence stars it is most likely not related to 
HESS Jl 626— 490. Furthermore, the upper limit for diffuse 
X-ray emission is a factor of ~2 lower than the expected 
value for HESS Jl 626—490 in a purely leptonic model. 

Not detecting any X-ray source fulfilling the energetic re- 
quirements for a purely leptonic scenario favors a hadronic 
emission process such as dense clouds in the vicinity of a 
cosmic particle accelerator. Such a scenario will be dis- 
cussed in the remaining part of this section. Dense molecu- 
lar clouds are established VHE 7-ray emitters because they 
provide target material in regions of high cosmic-ray den- 
sities. Using data from the Nanten 12 CO(J=1-0) Galactic 
plane survey, we detected a molecular cloud that is mor- 
phologically consistent with HESS J1626-490. This ob- 
ject is located at a kinematic distance of ~1.8 kpc. Fol- 
lowing [|T2l (Eqs. 2 and 3), we estimated the required 
gas density to produce the observed VHE 7-ray signal 
(F 7 (>0.6TeV) = 7.5 x 10- 12 ph cm" 2 s" 1 and Y = 2.2) as 
n « 10 cm -3 assuming a cosmic-ray production efficiency 
of 6 = 0.1 and a distance of d = 1.8 kpc. This value is an 
order of magnitude lower than the measured 12 CO mean 
density. Thus, this environment would be easily suited to 
providing the observed VHE 7-ray flux. 

Now assuming that this 12 CO / HI cloud is indeed the 
source of the observed VHE 7-rays, a cosmic-ray accel- 
erator would be needed in its vicinity. The nearby density 
depression seen in H I (see Sect. [3]) might indicate the pres- 
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Figure 4: Spectral energy distribution of HESS J1626-490 
showing the VHE 7-ray spectrum measured with H.E.S.S. 
(red). The solid line (blue) shows the hadronic model with 
the parameters that yield the best fit to the data (R = 35 pc, 

t = 1.5 x 10 5 yr, D 10 = 9 x 10 25 cm 2 s" 1 ). 

ence of a recent catastrophic event, such as an SNR, giving 
rise to strong shocks that would have blown the neutral gas 
out. At d — 1.8 kpc the edge of this region would be at a 
distance of 8.1 pc from the 12 CO/HI cloud. It is striking 
that this H I feature is consistent in both position and angu- 
lar extension with the SNRG335.2+00.1 0. Using the E- 
D relation the distance to SNR G335. 2+00.1 was estimated 
as d = 5 kpc |[T3lK which conflicts with the identification 
of SNR G335.2+00.1 with the H I depression seen at a dis- 
tance of 1.8 kpc. However, the scatter in S-D distances 
for individual sources is quite large (see [14']) making a 
positive identification of SNR G335.2+00. 1 and the H I de- 
pression seen at a kinematic distance of 1.8 kpc very well 
possible. 

To explore the physical scenario of the interaction of the 
SNR G335.2+00.1 with the 12 CO molecular cloud in more 
detail, we adopted the model described by |[T5l . In this 
model hadronic cosmic rays with a powerlaw energy spec- 
trum with index —2 and a total energy of 10 50 erg are 
injected at the SNR shock into the ISM. The transport of 
these cosmic rays is described by an energy-dependent dif- 
fusion coefficient of the form Dism = D10 • E - 5 , where 
E is the energy of the cosmic rays and Dio is the diffusion 
coefficient for an energy of 10 GeV. For the assumptions 
made by |[T5lL the transport equation can be solved analyt- 
ically and the expected 7-ray spectrum due to the produc- 
tion and subsequent decay of neutral pions can be calcu- 
lated. Some of the model parameters can be constrained 
by the measurements presented in this work. These are 
the mass M, the density n of the molecular cloud and the 
distance d between the object and the observer. Also, a 
lower limit for the physical distance R between the SNR 
shock and the cloud was derived, which corresponds to 
the projected distance between these two objects (8.1 pc). 
However, it is still possible that the actual separation be- 
tween SNR and cloud is in fact larger than this value due 



to slightly different line-of- sight distances between these 
two objects and the observer. The remaining free model 
parameters are the age of the SNR t, the diffusion coeffi- 
cient Dio and the radius R. Figure |4] shows the VHE 7-ray 
spectrum of HESS J1626— 490 measured with H.E.S.S., to- 
gether with the model with the parameters that yielded 
the best fit to the data (R = 35 pc, t = 1.5 x 10 5 yr, 
D 10 = 9 x 10 25 cm 2 s" 1 ). The best-fit value of D 10 is 
about two orders of magnitude lower than the estimated 
mean value for the Galactic plane, however, it agrees well 
with comparable measurements of diffusion in the vicinity 
of other cosmic-ray accelerators (e.g., for the W28 region, 
see lfT6l ). It has to be noted here, however, that the age 
t and the diffusion coefficient Dio are strongly correlated, 
and as t is not known for SNR G3 3 5. 2+00. 1 Dio is not very 
well constrained. Overall though, the modeling shows that 
a hadronic scenario is indeed viable for HESS Jl 626—490 
with reasonable physical parameters. Another complica- 
tion of hadronic models is the diffusion of cosmic rays into 
the target material which strongly depends on the density 
profile of the molecular cloud. More detailed and better re- 
solved molecular line observations could provide the nec- 
essary data in the future. 
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